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Abstract 



The effects of correlations on the bulk properties of nuclei are investigated 
in large model spaces including up to 21 single-particle orbits. The evaluation 
of the single-particle Green function is made feasible by means of the BAGEL 
approximation. The spectral function for single-nucleon pick-up and removal is 
investigated for the nuclei ^^O and '^^Ca. Special attention is paid to the effects 
produced by correlations on the calculated ground state properties of closed 
shell nuclei. It is observed that correlations beyond the Brueckner Hartree 
Fock approximation tend to improve the results obtained using realistic nucleon 
nucleon interactions. 
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1 Introduction 

The nuclear shell-model, which describes the nucleus as a system of nucleons moving 
without correlations in a mean field, has been very successful in explaining a large 
number of basic nuclear phenomena at low energies. Nevertheless, it has become 
quite obvious that to derive nuclear properties from a realistic nucleon-nucleon (NN) 
interaction it is necessary to consider the correlations between the nucleons. This 
necessity arises from the fact that NN interactions, which are adjusted to describe 
NN scattering data, contain very strong short-range components. As a consequence, 
Hartree-Fock (HF) calculations employing such interactions predict bulk properties of 
nuclei which are far off the empirical data and, in fact, fail to account for the binding 
of nuclei. 

Several methods have been developed to deal with the strong short-range com- 
ponents of the NN interaction. One of these methods is to consider that the wave- 
functions contain short-range correlations, e.g. Jastrow correlation functions, which 
are optimized in a variational calculation An alternative approach is to introduce 
effective operators. A typical example of such an effective operator is the Brueckner 
G-matrix |Q, |^. This G-matrix is evaluated from the NN interaction by solving the 
Bethe-Goldstone equation, which accounts for the virtual excitation of two interacting 
nucleons into single-particle states above the Fermi energy, i.e states unoccupied by 
other nucleons. In solving the Bethe-Goldstone equation for realistic NN interactions 
one has to consider 2-particle configurations up to an excitation energy of a few GeV 

The consideration of the short range, or high-energy, correlations in the Brueckner- 
Hartree-Fock (BHF) approach results in a major improvement on the prediction of 
nuclear properties, as compared to the independent-particle or HF approach. How- 
ever, trying to evaluate the binding-energy and saturation density of nuclear matter, 
or the binding-energies and radii of finite nuclei within the BHF approximation, one 
still observes significant deviations from the corresponding data Motivated by 
the success of the Walecka model [Q , attempts have been made to include relativistic 
effects within the so-called Dirac-BHF approach, leading to quite an improvement for 
the description of nuclear matter 0. The corresponding calculation for finite nuclei, 
however, did not yield a satisfying agreement between theory and experimental data 

It is worth noting that quite different results are obtained in the studies of infinite 
nuclear matter and finite nuclei. This demonstrates that the structure of finite nuclei 
is more complicated and cannot be treated simply as a piece of nuclear matter of finite 
size. These differences cannot be attributed to the short-range correlations, as the 
configurations at high energies should be similar in both cases. It is, therefore, more 
likely that the different behaviour of nuclear matter and finite nuclei is due to the 
excitation modes of these systems at lower energies. At such energies the excitation 
spectrum for a finite system still reflects the shell structure, which is absent for the 
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infinite system. 

In this study, we investigate improvements to the BHF approach caused by a 
detailed treatment of low- and medium-energy excitations. Specifically, we investigate 
the effects produced by configurations with excitation energy up to about 100 MeV. 
These configurations are formed from single-particle states which all belong to a few 
major shells around the Fermi level. The effects produced on the nuclear wavefunction 
by the admixtures of these configurations will be described in the following as long- 
range correlations, in contrast to the short-range correlations discussed above. 

For the study of these long-range correlations we employ a Green-function ap- 
proach within a finite model space of discrete single-particle states [|, |TT], [T^, |T3|. We 



describe the irreducible self-energy for the nucleons in terms of the HF contribution 
plus terms of higher order in the residual NN interaction, which are described by the 
diagrams shown in Fig. 1. For the residual interaction Q we consider a G- matrix, 
which sums up all contributions originating from high-energy 2-particle excitations, 
i.e excitations outside the space of our calculation (model-space). Thus all intermedi- 
ate states displayed in Fig. 1 must be understood to belong to the model-space under 
consideration. Our definition of the nucleon self-energy is different from that em- 
ployed in the BHF approximation, which includes diagram la, and those containing 
particle-particle ladder diagrams like Ic, but ignores contributions like those described 
by diagrams lb. Id and le. This non-symmetric treatment of particle-particle and 
hole-hole excitation is a typical feature of the Brueckner hole-line expansion. It may 
very well be justified for the study of short-range correlations, as the high-energy 
particle-particle configurations do not have any counterpart on the hole-hole side. 
The ordering of the hole-line expansion, however, may be questionable if low-energy 
excitations are to be considered. Therefore we would like to consider in our approxi- 
mation for the self-energy, the contribution of all terms with intermediate 2 particle-1 
hole and 2 hole-1 particle states (see Fig. 1) including all interactions between those 
states. 

For the irreducible self-energy one may then solve the Dyson equation to deter- 
mine the reducible self-energy or the corresponding single-particle Green-function. In 
the Lehmann representation this Green function is defined in terms of many poles 
below and above the Fermi-energy and the value of the spectral function at these 
poles. In ref. fl^ the so-called BAGEL approach has been introduced to represent 
the complicated structure of the Green-function in terms of a few "characteristic" 
poles. Using this approximate representation, it becomes possible to evaluate a self- 
consistent Green-function in the sense that the irreducible self-energy is calculated 
using the Green-function, which, in turn, is obtained from the Dyson equation for 
this self-energy ||14 |. 

One of the aims of the present study is to investigate the sensitivity of the results 
for the Green-function or the spectral function on the approximations chosen for the 
irreducible self-energy and on the choice of the model-space. As examples we consider 
the closed shell nuclei ^^Oand ^°Ca. The results obtained for the spectral function 
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can directly be compared with experimental data obtained in (e, e')p experiments 



[p!5| , |T7|. From the Green- function we can furthermore derive the occupation 
probabilities for the various single-particle orbits. These occupation probabilities are 
quite a conspicuous measure for the deviation from the independent particle model, 
i.e. a measure for the amount of correlation effects. Such occupation probabilities 
may be compared to results obtained for nuclear matter by using the Green-function 
approach |18[, within Brueckner theory [19, 20 1 or by means of the correlated basis 
function theory ||2l |. 

Finally, one may also use the Green-function to determine the ground-state energy 
and other observables like the radius. One can study the effects of correlations beyond 
the BHF approximation and estimate the importance of particle-particle and hole- 
hole ladder diagrams on the saturation properties of infinite nuclear matter [^, 

The Green-function approach is reviewed briefly in sect. 2. The presentation there 
helps to establish the nomenclature but also to define some details of the present 
investigation. Results on the spectral function, occupation probabilities, binding 
energies and radii for ^^O and ^^Ca obtained with the use of various approximations, 
are presented in sect. 3. The main conclusions are summarized in sect. 4. 



2 BAGEL approximation for the one-body Green- 
function 

2.1 Model Space and effective Hamiltonian 

As already mentioned in the introduction, one of the main aims of the present investi- 
gation is to explore effects of long-range correlations by means of the Green-function 
approach within a certain truncated model-space. On the other hand, the effects of 
short-range correlations are taken into account by the introduction of an effective 
interaction, i.e a G-matrix appropriate for the model-space. Such a concept of a 



double-partitioned Hilbert space has been used before for finite nuclei [O, El] and 



infinite nuclear matter [^, The G-matrix is determined as the solution of the 



Bethe-Goldstone equation 

G = v + v — g. (1) 

^ "^mod^ "^mod 

In this equation T is the kinetic energy operator while V stands for the bare two- 
body interaction. For the latter we have chosen the One-Boson-Exchange potential 
C defined in table A.l of The Pauli operator Q^aQd is defined in terms of 



harmonic-oscillator single-particle states. Thus applying Qj^^od two-particle states 
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I a/3 > one obtains 

{0 if a or /3 below Fermi level 

if a and (3 in model space (2) 

I a/5 > else 

Here the Fermi level refers to the Fermi level of the independent particle model for the 
nucleus under consideration. The model-space in eq. consists of all single-particle 
states up to and including the 0hlf2p shell of the oscillator potential, while a constant 
value of = —30 MeV has been adopted for the starting energy. This choice of using 
a constant starting energy in the G-matrix is of course an approximation introduced 
to simplify the calculations. In sect. 3 we discuss the sensitivity of our results on the 
choice of the starting energy. There we will also consider the effects of using different 
model spaces. 

A major problem in the use of the modern OBE versions of realistic NN inter- 



actions, which are defined e.g. in [^, is the fact that BHF calculations using these 
interactions yield predictions for the radii which are too small compared to the ex- 
perimental data 1^ . This has serious consequences on calculations beyond the mean 
field approximation. If one uses e.g. such BHF wave functions to evaluate the resid- 
ual interaction to be used in shell-model calculations, one obtains results which are 
in poor agreement with the experimental data. On the other hand, evaluating the 
matrix elements for a set of harmonic oscillator wave functions which yield accurate 
radii, leads to matrix elements for the shell-model calculation, which provide a fair 
agreement with experiment [^. Therefore one is tempted to modify the hamiltonian 
in such a way that a HF calculation using Q yields a self-consistent basis of oscillator 
wave functions which reproduce the experimental radius. This is achieved by sub- 
tracting from the Q defined in eq. (|1|) the non-diagonal terms of the HF single-particle 
hamiltonian in the oscillator basis 

= ^ {T^p + Uafi) alttfs , (3) 

where Tq,^ and Uais define the matrix elements of the kinetic energy and the HF 
single-particle potential , respectively, calculated in a harmonic oscillator basis. This 
basis is defined by an oscillator parameter of 6 = 1.76 fm for all studies of ^^O and 
b = 2.0 fm in case of the nucleus ^°Ca. 



2.2 Green function and matrix diagonalisation 

As a first step towards the evaluation of the single-particle Green function we deter- 
mine the contributions to the irreducible self-energy, which are of second order in the 
residual interaction: 



(4) 



■ySfj. ' 
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Note that the summation on single-particle states 7 , 5 , /i is restricted to states within 
the model-space. For discrete states a, the ga{^^) denote the Hartree-Fock (HF) 
approximation for the Green function |T3]. This equation can be rewritten into the 
form 

'^Pi^2,h uj-eipi,p2,h)+t7] '^h5^,,p uj-e{hi,h2,p)-iv 



(5) 



where we have introduced the abbreviation 



e(a,/3,7) = 6r + er-<'^ (6) 

and the e^^ are the HF single-particle energies. In eq. (^) the summations on particle 
labels like pi,P2 and p are restricted to those single-particle states within the model 
space, which are above the Fermi level, whereas the labels hi, /12 and h refer to hole 
states. It is evident that the first term on the right hand side of eq. @ refers to the 2 
particle - 1 hole contribution to the self-energy (diagram a of Fig. 1) while the second 
term defines 2 hole - 1 particle contribution (diagram lb). This self-energy can now 
be inserted into a Dyson equation for the Green function Gai3{uj) taking into account 
the correlation effects contained in E^^^ 

Ga/si^) = Sal3ga{uj) +J2-y9a{i^)^a'Y{i^)G^f3{uj) (7) 



^7 • 



Lj—u)„ +iri 



The second line (eq.(|])) exhibits the Lehmann representation of the single-particle 
Green function in terms of the spectroscopic amplitudes < \l/^|aQ|\l/^+^ > and 
< \l/^+^|a)j|\l/o >, where Qa (0)3) stands for the single-particle annihilation (creation) 
operator in the HF basis. The state refers to the ground-state of the A-particle 
system, while ^l/^'*'^ (^m~^) stands for the states of the A+l-particle (A-1) system as 
obtained in the present approach. Also the energy variables 

=E^^'-E^. (9) 

= Eq - E^~^ , (10) 



are defined in terms of energies obtained for the states of the nuclei with A, A + 1 
and A — 1 nucleons. The substantial ingredients of the Lehmann representation of 
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eq. (H) can be obtained from a solution of an eigenvalue problem |T3[ 
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(11) 



In writing this equation, we assume that for a given set of conserved quantum numbers 
(parity, isospin, angular momentum) we have in our model space two HF single- 
particle states (a, P) and K 2plh configurations with energies (see definition of 
eq. i)) 

ei = e{pi,p2,h) , (12) 



and connecting matrix elements 



b,. 



=< ah\g\pip2 > 

--< ms\PiP2 > 



(13) 



Furthermore we consider L 2hlp configurations of the same symmetry as the single- 
particle states. These have energies 



and connecting matrix elements 

A, 

B., 



Ej = e{hi, h2,p) 



=< ap\Q\hih2 > 

--<pp\g\hih2 > 



(14) 



(15) 



Solving eq. ([11]) one obtains eigenvalues Un- These eigenvalues must be identified 
with uj^ of the Lehmann representation of eq. (|^), if Un is an energy above the Fermi 
energyQ In this case we have 



If Un is below the Fermi energy we obtain Un = uj~ and 



(16) 



(17) 



In a straightforward way one can improve the approximation discussed so far and 
incorporate the effects of residual interactions between the 2plh and 2hlp configu- 
rations, i.e. account for diagrams like those displayed in Fig. Ic, Id and le in the 

"'^The definition of the Fermi energy is given below 
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definition of the self-energy. One simply has to modify the corresponding parts of the 
matrix in eq. ([TT|) and replace 

/ei ... \ /El ... \ 

: ■ • . \ ^ ^2pih , and i ' • . \ ^ ^2^1^ , (18) 

\0 ... ck) VO ... Ek) 

where TC2pih and H2hip contain the residual interactions in the 2plh and 2hlp sub- 
spaces. 

The normalization of the eigenvectors of eq. ( [TT| ) ensures that 



n m 



where the sum on n accounts for all solutions with larger than the Fermi energy 
and the sum on m for all solutions with energies a;~ below the Fermi energy. These 
spectroscopic factors can also be used to determine the occupation probabilities 



N^ = J2\<^i-'K\^o>\' ■ (20) 

m 

The approach discussed up to now, however, is not number-conserving in the sense 
that it is not guaranteed that 

EiVa(2ja + l)=A, (21) 

a 

with {2ja + 1) representing the degeneracy of the states a. This equation, however, 
allows one to determine the Fermi energy in such a way that the particle number 
calculated according to the left hand side of eq. ( pT]) is as close as possible to the 
mass number A of the nucleus under consideration. An approach which is strictly 
number-conserving is obtained only if the Green function used to calculate the self- 
energy (like in eq. (^) is consistent with the Green function resulting from the solution 
of the Dyson equation (0) 0. 

With the help of the single-particle Green function one can evaluate the 
expectation value for an arbitrary single-particle operator O 



< ^o|0|^o >= / ^ E < > G^piu) , (22) 

Jc 2m 



where the C below the integral sign denotes a contour integration counter-clockwise 
in the upper half plane including the real axis. Therefore only the contributions from 
the poles at energies below the Fermi energy have to be considered and using the 
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nomenclature of the matrix representation in eq. (0) this expectation value can be 
rewritten as 

< ^o\0\^o >=J2<^\0\P> X^^X^p , (23) 

af3m 

where the sum on m is restricted to solutions with energies a;" below the Fermi 
energy. The Green function can furthermore be used to evaluate the binding energy 
of the A-nucleon system by 




where T denotes the kinetic energy. 



2.3 The BAGEL approximation 



After we have rewritten the solution of the Dyson equation in terms of an eigenvalue 
problem (eq. (|TT|) ), it is now straightforward to introduce the BAGEL approximation. 



The BAGEL approximation has originally been formulated to define an appropri- 
ate truncation scheme for shell-model calculations in large model spaces. Since then 



this method has been employed in several calculations |31], including attempts to 



derive effective operators for truncated model spaces |^2[. The manner in which the 
BAGEL approximation can be used to determine the single-particle Green function 



has already been discussed in [0, , so that only a brief outline of the method need 
to be given here. 

For many applications of the single-particle Green function it is not necessary to 
know in detail all of its poles in the Lehmann representation. It should rather be 
sufficient to describe the energy distribution of the spectroscopic amplitudes in terms 
of a few "characteristic" states of the A ± 1 systems. For that purpose we consider 
the operator a which corresponds to a part of the matrix in eqs. (|rT|) and 



ai 



H 



2plh 



J 



(25) 



and apply this operator on the single-particle state \a >, which in terms of the matrix 
representation of eq. (|25|) is described by the column vector (1, ... 0)^ 



a\a >- 



e \a > +ai «! > 



(26) 



where \ai > is orthogonal to |a > and the coefficient di is chosen so that \ai > is 
normalized. Following the Lanczos algorithm one can subsequently construct 
additional states >, which are all orthogonal to each other. Applying the Lanczos 
procedure N times one obtains N basis states of the 2plh configuration space, which 
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have the same symmetry quantum numbers as the single-particle state \a >. If now 
the model-space contains another single-particle state |/3 > of the same symmetry, 
one may consider in an analogous way the operator b defined by 



bi 



n 



2plh 



V b 



K 



bK\ 
J 



(27) 



and obtain by applying it initially to |/3 > additional N basis states \Pi >■ In this 
case special care should be taken that the states |/3j > are not orthonormalized only 
among themselves but also with respect to the \ai >. 

In a similar way we can furthermore construct for each single-particle state \a >, 
M basis states of the 2hlp configuration space by considering the corresponding sub 
matrices of eg. (pH]) 

Ai ... Al\ 



HF 



A 



Ai 



\Ar 



n 



2hlp 



{2i 



) 



In this way, the BAGEL(N,M) approximation reduces the matrix defined in eq.([TTD 
to a subspace of dimension / * [N + M + 1), if / is the number of single-particle states 
of the same symmetry. The Green function of this BAGEL (N,M) approximation is 
then defined in the manner discussed in the preceding section but now considering 
only the eigenvalues and eigenvectors of the truncated matrix. It is obvious that the 
BAGEL (0,0) corresponds to the HF approximation, while for large enough N and M 
the BAGEL (N,M) becomes identical to the exact solution of eq. (|TT|). 



3 Results and Discussion 
3.1 The Spectral Function 

The spectral function summarizes the spectroscopic factors for adding, at energies 
above the Fermi energy Ep, or removing, at energies below Ep, nucleons from the 
nucleus under consideration. For a model space with discrete single-particle orbits, 
as we consider in the present investigation, the spectral function is given in terms of 
5 functions at the various poles of the single-particle Green function (see eq. (|D). For 
the sake of presentation, we will fold these 5 functions with a gaussian distribution 
assuming a with F. This width can be interpreted as a simple way to include the 
escape width for the states of the residual nuclei, or as a tool to express the uncertainty 
of the theoretical calculation. The second interpretation applies in particular to the 
use of the BAGEL approximation. With this folding procedure the spectral function 
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is given by 



+ E„^exp-(^^^^| < ^t'Kl'^o > ?] . (29) 



In this equation, as in eq. (|I^), the summation over n is restricted to the poles 
above Ep and the summation over m to those below the Fermi energy. 

As a first example, we present in Fig. 2 results for the spectral function of the 
P3/2 orbits in '^^Ca for energies below the Fermi energy. In this figure we compare 
the results obtained in the BAGEL(M=10, N=10) approach to those obtained using 
M=N=20. If we choose for the presentation a width of F = 2 MeV (lower part of 
the figure) differences between the two approximations are clearly visible. It is not 
only that the M=N=20 distribution displays more poles, as the corresponding Green 
function contains more poles, but also the position of the peaks are shifted if the 
accuracy of the BAGEL approximation is improved. If, however, we examine the 
same results with a reduced energy resolution, obtained by increasing the width in 
eq.(P^) to F=5 MeV, the differences between the 2 BAGEL approximation become 
invisible. This is a demonstration of the fact that the BAGEL(M,M) approximation 
reproduces all energy-weighted moments of the exact calculated distribution from 



order n = to the order n = 2M + 1 |T4| . In the following discussion, we will present 
results for the BAGEL (20,20) approximation and display the data assuming a width 
F=2MeV. 

The spectral functions displayed in figures 3 and 4 have been obtained by mul- 
tiplying the function Sa above with the degeneracy factors 2ja + 1 and adding up 
the various contributions with the same orbital angular momentum. Fig. 3 shows 
the spectral functions for ^^O in the energy interval ranging from uj = -50 MeV to 
30 MeV. The results in this figure clearly reflect the shell-structure of nuclei and the 
limitations of the independent particle model (IPM). For orbital angular momentum 
/ = one observes a strong peak at -4 MeV, just above the Fermi energy, which can 
be identified with the lsi/2 state of the IPM. The single-particle strength for the Osi/2 
and 2si/2 states is distributed in energy intervals around -40 MeV and around +20 
MeV, respectively. 

Similar features are also observed for the / = 1 states. The strength below the 
Fermi energy is mainly located at an energy u around -18 MeV. This peak originates 
from two poles: one at a;=-19.12 MeV (P3/2), the other at -16.73 MeV (^1/2)- Each of 
these poles carries around 75% of the total strength below the Fermi energy in P3/2 
and pi/2 states, respectively. With the energy resolution (F = 2 MeV) considered 
for the figures, the spin-orbit splitting is not resolved. It is worth noting that the 
contributions to the self-energy beyond HF reduce the calculated spin-orbit splitting 
for the / = 1 hole states from 3.6 MeV to 2.4 MeV. This last number is calculated 
by taking the energy difference between the two states carrying maximum strength. 
The decrease in the spin-orbit splitting is in agreement with the results of previous 
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calculations |^4[, which suggest that relativistic effects are mainly responsible for the 
spin-orbit splitting between hole states. 

Additional / = 1 single-particle strength is observed at lower energies, although the 
individual contributions are so small that they are barely visible in Fig. 3. It is worth 
noting that the total occupation probability (see eq. (pOD) is 0.926 (0.901) for the P3/2 
(Pi/2) state while the contribution of the main peak to this value is only 0.781 (0.790). 
The additional strength originates from states at lower energies, which implies larger 
excitation energies in the A-1 nucleus. The / = 1 spectral function furthermore 
displays the single-particle strength at positive energies, which corresponds to the Ip 
and 2p states of the IPM. 

In the top part of Fig. 3 the spectral function for the / = 2 and / = 3 states are 
displayed. In the IPM these states are unoccupied. Due to the correlations included 
in our approximation, we observe an occupation probability of 0.052 and 0.010 for the 
d and /-states, respectively. In the case of the / = 2 spectral function one can observe 
two peaks originating from the ^5/2 and ^3/2 states. The calculated spin orbit splitting 
between these states is 4.92 MeV, slightly below the HF result of 5.41 MeV. For these 
particle states, the result with and without correlations, is in good agreement with 
the experimental value of 5.09 MeV. This is again in line with the arguments given in 
ref. IQ, as the relativistic effects discussed there should not influence the spin-orbit 
splitting for the particle states. 

In Fig. 4 we show the spectral function calculated for ^°Ca, considering the orbital 
angular momenta / = to / = 3. The presentation is limited to energies below the 
Fermi energy, since for these energies a comparison with experimental data, deduced 



from (e,e')p experiments is possible |T7]. In the case of / = 0, the spectral strength for 
the hole-state close to the Fermi energy is nicely reproduced. The broad distribution 
of / = strength originating from deep-lying hole states, which is predicted in the 
calculation down to energies of -100 MeV, has not been resolved in the experimental 
data. Also for ^°Ca the spin-orbit sphtting calculated for the d hole states is smaller 
than the observed splitting. In the experimental data one finds a spectral strength 
below the Fermi-energy for / = 3 (notice the modified scale in the lower right part of 
Fig. 4). Such a strength is completely absent in the IPM or HF approach. Including 
correlation effects we observe a strength of similar size as the experimental data. 
While, however, the experimental strength is located in essentially one state, the 
calculation yields a distribution over various states. The agreement between theory 
and experiment is reasonable, keeping in mind that the theoretical results have been 
derived from a realistic NN interaction without any adjustable parameter. Results 
on the occupation probabilities obtained for ^^Ca considering various approximations 
for the nucleon self-energy are listed in table |l[ 
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3.2 Correlations and Ground-State Properties 

A major aim of the present investigation is to study the effects of correlations included 
in the definition of the nucleon self-energy on the bulk properties of closed shell nuclei. 
For that purpose we consider various approximations. The first of these is what we will 
call the Hartree-Fock (HF) approximation, i.e. the self-energy is calculated in the HF 
approximation. Note, however, that this HF approximation is something in between 
a conventional HF and BHF calculation. The G-matrix interaction (see section 2.1) 
accounts for intermediate 2-particle states, which are outside the considered model 
space. However, in contrast to a self-consistent BHF calculation it does not account 
for 2-particle states between the Fermi level and the limit of the model space. Such 
correlations are taken into account in the subsequent approximations. Furthermore 
the hamiltonian has been adjusted in such a way that this HF approximation leads to 
a value for the radius of the nuclei ^^O and ^°Ca, which are close to the experimental 
values (eq.(ID). Therefore, in studying effects of correlations beyond HF, one should 
pay attention only to the effects in the calculated radii relative to the HF result and 
not expect any improvement compared to the experiment. 

The second approximation considered in our study consists in adding the con- 
tribution of the second-order (2plh) diagram (a) in Fig. 1 to the self-energy. This 
approximation will be labeled "bare 2plh" approach in the tables and the discussion 
below. A part of the residual interaction in the space of 2plh configurations is given 
by the interaction between the two particle lines. The approximation which accounts 
for this part of the residual interaction and includes the diagram of Fig. Ic plus all 
other particle-particle ladder diagrams in the definition of the self-energy, corresponds 
to the Brueckner - Hartree-Fock approach and will be identified in the following as 
"BHF" . On the other hand the term "2plh" describes the approximation in which the 
whole residual interaction is considered in the 2plh subspace (diagram of Fig. Id and 
others). Finally, the term "total" describes the calculation in the complete space of 
2plh and 2hlp configurations with all residual interactions included. This means that 
the complete matrix in eq.(|ll]) extended by the modifications of eq.(0) is considered 
to construct the Green function. 

For all studies beyond HF the single-particle Green function has been determined 
in the BAGEL(M=10,N=10) approximation. The results on the bulk properties of 
nuclei discussed below turned out to be extremely stable with respect to this approx- 
imation. Increasing the accuracy of the BAGEL to M=N=20, did not lead to any 
change in the results displayed in the tables. 

For the resulting Green functions we have determined the occupation probabilities 
Na according to eq. (l20| ) and a mean value for the spectral distribution below the 
Fermi energy, defined by 

The contour integral in the first part of the equation is defined in the same way as in 
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eq. p^ ) and as before the summation in the second part of the equation is restricted to 
poles in the Green function with an energy u;" below the Fermi energy. In the case of 
the HF approximation the definition in eq. ( PPD leads to the HF single-particle energy. 
The modification of these mean values due to the correlations give some insight into 
the effects of correlation on the calculated binding energy (compare eq.(p^). 

Results for the occupation probabilities and the e< are listed in table |I] for some 
single-particle orbits of ^°Ca. One finds that the admixture of 2plh configurations 
reduces the occupation probability by 3% for the deep-lying hole states and by around 
5% for states closer to the Fermi level. Inclusion of the residual interaction between 
the 2plh configurations slightly enhances these depletions for the occupation probabil- 
ities. The effects of the residual interaction are dominated by the pp ladder diagrams 
while the additional interaction terms can essentially be ignored. The removal of oc- 
cupation strength from the state below the Fermi energy to 2plh configurations above 
Ep is accompanied by a lowering of the mean values e< below the HF single-particle 



e . This is obvious since one has the following relation [ffS 



Y: iVae<,a + (1 - iVa)e>,a = E ^ ' (31) 
a a 

where the sum on a is a sum on all single-particle orbits of the same symmetry and 
the mean value for the spectral distribution above the Fermi energy is defined in 
analogy to eq. (|30D by 

' ■^^n+i<^r^i4i^o>r- (32) 



It is this lowering of the mean values e<, which is the main source for the gain in 
energy per particle when we evaluate the energy according to eq. (|2^ ) using the "bare 
2plh", "BHF" and "2plh" approximations. The calculated values for the energy per 
particle for ^^O and ^^Ca are listed in table ^ As seen in this table, the inclusion of 
2plh configurations within our model space increases the calculated binding energies 
by 1.9 MeV and 1.6 MeV per nucleon for ^^O and ^°Ca, respectively. As already dis- 
cussed above, the effects of the residual interaction between 2plh configurations, are 
small (0.3 MeV per nucleon) and are mainly due to the particle-particle interactions 
contained in the BHF approach. The effect of the 2plh admixture on the calculated 
radii is essentially negligible. 

It should be mentioned that the way we calculate the energy per nucleon in the 
BHF approximation is slightly different from the conventional form. In conventional 
BHF calculations one determines the total energy from the mean values e<, which, for 
the corresponding approximation in the self-energy, are identical to the BHF, e^^^ , 
single-particle energies and calculates the total energy from 



pBHF 

conv 



^E[<«l^l">+ef''^]^"' (33) 
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with Ua the occupation probabihties of the HF or IPM approach. This expression 
essentially corresponds to eq. (24) except that the conventional BHF approach ignores 
the differences in the occupation probabilities for the various states. The calculated 
binding energy per nucleon according to eq. (p^), using for the particle number the 
results of eq. (0) are very close to those obtained in the conventional way (differences 
are typically around 0.1 MeV). 

The inclusion of 2hlp configurations in the "total" calculation leads to a frag- 
mentation of the strength into many states below the Fermi energy. One finds, that 
occupation probabilities and mean values for the spectral distribution e< are not 
very much affected by the inclusion of these additional correlations (see table p for 
the states which are occupied in the IPM. The total approach, however, leads to a 
non- vanishing occupation also for configurations with quantum numbers a, which are 
unoccupied in the IPM. The mean value e< for these 2hlp configurations is typically 
quite attractive. Therefore one observes a non-negligible gain in the binding energy 
due to the hole- hole correlations. In our standard model space this gain is 0.4 MeV 
and 0.3 MeV per nucleon for ^^O and '^^Ca, respectively. 

Considering just the energies it is obvious, that the effect of 2plh correlations 
is significantly larger than those obtained from 2hlp. This would support the BHF 
approximation, which includes the effects of the former kind but ignores 2hlp contri- 
butions. For the calculation of the nuclear radii, however, the 2hlp correlations seem 
to be more important than those deduced from particle-particle correlations. The 
inclusion of 2hlp correlations yields an increase for the calculated radii of around 0.1 
fm. Alone this increase does not seem to be dramatic. It is worth noting, however, 
that the combination, an increase of calculated binding energy and radius, moves the 



result for these bulk properties of nuclei away from the so-called "Coester band" p7 |. 
This is a similar feature as has also been observed for hole-hole correlations in infinite 
nuclear matter [^, 

The sensitivity of our results on the choice of the model space and on the param- 
eters of the residual interaction is examined in table |^. This table contains results 
for ^^O for three choices of the model space: one with 10 active orbits, another with 
15 and a third where 21 active orbits are considered. This table contains results for 
the energies per nucleon obtained in the HF approximation, for the "total" approach 
and the difference between these 2 values. Since the effective interaction depends on 
the choice of the model space, also the HF results change significantly. For the small 
model space {Ncut=^0) the effects of 2plh correlation with 2 particles in sdg-shell 
or above are contained already in the interaction. Therefore a larger binding energy 
is obtained in the HF approximation for this small model space than for the HF 
calculation using the interaction which is appropriate for a larger model space. 

These effects of the 2plh correlations and the interference with 2hlp correlations 
are contained in the "total" calculation. The convergence of these correlation effects 
is rather slow with increasing model space. Both the calculated binding energies as 
well as the results for the radii are enlarged if the number of shells included in the 
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model space is increased. Therefore the features discussed above may get enhanced 
if even larger model spaces are considered. 

Table ^ displays results also for an alternative choice for the starting energy in 
eq. (|l]), using u = —5MeV instead of our standard choice u = —30MeV. One 
observes that the larger value for u yields larger energies already in the HF approx- 
imation and also for the "total" calculation. The correlation effects reflected in the 
energy differences between these two approximations are rather insensitive to the 
choice of u. 



4 Conclusions 

The effects of correlations beyond the Hartree-Fock approximation are investigated 
in large model spaces including up to 21 active single-particle states of the harmonic 
oscillator potential. The residual interaction for this model space is derived from a 
realistic meson-exchange potential accounting for short-range correlation effects 



beyond the model space by solving an appropriate Bethe-Goldstone equation. The 
correlations within the model space are evaluated in terms of the single-particle Green 
function. The Green function can easily be evaluated even for a large model space 



by employing the BAGEL [[131 , |T4| approximation. 

The spectral functions evaluated from this Green function for ^^O and ^^Ca yield a 
distribution of single-particle strength over a broad range of energies. A comparison 
with an empirical spectral function for "^^Ca, determined from (e, e')p experiments 
1^ demonstrates that the spectral function, calculated without adjustable parame- 



ter, produces a reasonable agreement. However, the spin-orbit splitting between hole 
states close to the Fermi energy is underestimated even after the inclusion of correla- 
tion effects. It is argued that relativistic features, not considered in the present work, 
are mainly responsible for this splitting p| . 



Special attention is paid to the effects of correlation, treated in various approxima- 
tions, on the calculation of bulk properties of nuclei, like binding energy and radius. 
It turns out that the effects of correlation on the calculated energies are dominated by 
the particle-particle (2plh) correlations taken into account in the BHF approxima- 
tion. Some additional binding energy is obtained if also hole-hole (2hlp) correlations 
are considered. The 2hlp correlations furthermore give rise to an increase for the 
calculated radius. Therefore the properties evaluated for ground states of nuclei are 



"moved off the Coester band" towards the experimental data |2^. This effect of 
hole-hole correlations is similar to the one observed in nuclear matter ^ . 

The effects obtained from inclusion of the hole-hole correlations are not very large. 
In the model space under investigation they yield an extra gain in binding energy of 
around 0.4 MeV per nucleon and an increase of the radii of around 0.1 fm. The effects 
may get larger, if larger model spaces are taken into account. This is possible with 
the techniques explored here. However, for a realistic investigation one should try to 
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represent the particle states in terms of plane waves rather than harmonic oscillator 
states. 

This work has partly been supported by the Deutsche Forschungsgemeinschaft 
(DFG, Fa67/14-1). 
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Table 1: Occupation probabilities and mean values for the energy distribution below 
the Fermi energy, calculated for different single-particle states in ^^Ca. The differ- 
ent columns correspond to various approximations used in the calculation for the 
self-energy of the nucleons. The column labeled "HF" stands for the Hartree-Fock 
approximation, whilein the second column also the second-order diagram of Fig. la has 
been taken into account ("bare 2plh"). All pp ladder diagrams (like e.g. the diagram 
of Fig.lc) are included in the Brueckner-HF ("BHF") approximation. If all residual 
interactions between 2plh configurations (also Fig. Id) are taken into account, one 
obtains column 4 ("2plh"). The results of the complete calculation, also including 
the 2hlp configurations with residual interactions are given in column 5 ("total"). 
The single-particle energies are given in MeV. 





HF 


bare 2plh 


BHF 


2plh 


total 


Osi/2, e< 


-58.211 
1.0000 


-61.098 
0.9722 


-61.189 
0.9704 


-61.126 
0.9714 


-61.368 
0.9689 


0p3/2, e< 


-38.796 
1.0000 


-41.552 
0.9672 


-41.695 
0.9637 


-41.649 
0.9642 


-41.876 
0.9610 


Opi/2, e< 


-35.901 
1.0000 


-39.074 
0.9618 


-39.236 
0.9577 


-39.165 
0.9589 


-39.432 
0.9553 


04/2, e< 


-20.747 
1.0000 


-23.597 
0.9511 


-23.905 
0.9383 


-23.941 
0.9335 


-24.195 
0.9267 


04/2, e< 

Na 


-19.446 
1.0000 


-19.861 
0.9363 


-20.231 
0.9206 


-20.170 
0.9199 


-20.509 
0.9112 


lsi/2, e< 

Na 


-16.350 
1.0000 


-22.661 
0.9439 


-23.003 
0.9298 


-23.009 
0.9259 


-23.281 
0.9187 


0/7/2, e< 

Na 


0.0 


0.0 


0.0 


0.0 


-40.058 
0.0459 


0/5/2, e< 

Na 


0.0 


0.0 


0.0 


0.0 


-42.952 
0.0444 


lP3/2, e< 

Na 


0.0 


0.0 


0.0 


0.0 


-48.003 
0.0381 


lPi/2, e< 

Na 


0.0 


0.0 


0.0 


0.0 


-48.254 
0.0391 
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Table 2: Results for binding energy per nucleon (E/A), radius of the nucleon distribu- 
tion and particle number A obtained in various approximations (see caption of table 
1). The effects of the Coulomb energy and the charge radius of the proton have been 
removed from the experimental data to allow a direct comparison to the calculated 
values. Results are given for ^^O and "^^Ca. For the calculation of ^^O particle states 
up to and including the sdg shell arc taken into account, while in the case of ^^Ca all 
states up to the hfp shell are considered. All energies are listed in MeV and the radii 
in fm. 





WQ 

E/A Radius A 


40Ca 

E/A Radius A 


HF 
bare 2plh 
BHF 
2plh 
total 


-4.202 2.633 16.000 
-6.099 2.626 15.159 
-6.373 2.625 14.894 
-6.370 2.625 14.868 
-6.747 2.720 16.262 


-6.491 3.455 40.000 
-8.104 3.451 38.152 
-8.275 3.449 37.765 
-8.278 3.448 37.700 
-8.587 3.533 40.379 


Exp 


-9.118 2.580 16.000 


-10.597 3.410 40.000 
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Table 3: Results for the energy per nucleon in the HF approximation, in the "total" 
approximation and the gain in energy AE relative to HF, obtained in various model 
spaces for the nucleus ^^O. Furthermore the radius of the nucleon distribution and 
the particle number A are listed calculated in the "total" approximation (see tables 
1 and 2). The various model spaces are identified in terms of Ncut the number of 
active single-particle orbits. 



Ncut 


Ehf/A 


Etot/A 


AE/A Radius 


A 








=-30 MeV 






10 


-4.791 


-6.388 


-1.597 


2.683 


16.155 


15 


-4.202 


-6.747 


-2.545 


2.720 


16.262 


21 


-3.628 


-7.052 


-3.424 


2.754 


16.379 






UJ 


=-5 MeV 






10 


-6.455 


-8.171 


-1.716 


2.686 


16.162 


15 


-5.556 


-8.214 


-2.658 


2.723 


16.279 


21 


-4.863 


-8.398 


-3.535 


2.756 


16.384 
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5 Figure Captions 

Figure 1: Diagrams representing various contributions to the self-energy of second 
and higher order in the residual interaction Q. 

Figure 2: The spectral function defined in eq. (29) obtained for the P3/2 states in 
^°Ca, employing two different BAGEL approximations. The curves in the upper part 
are obtained assuming a width of r=5 MeV, whereas in the presentation displayed 
in the lower part r=2 MeV has been used. 

Figure 3: Spectral functions for the nucleus in a region of energies between 
a;=-50 MeV to lo— 30 MeV. The spectral functions for the various orbital angular 
momenta {I — 0, 1, 2, 3) have been obtained from eq.(29) (width r=2 MeV) by mul- 
tiplying the contributions from the different total angular momenta = Zq, ± 1/2 by 
the degeneracy factors 2ja + 1. 

Figure 4: Spectral functions obtained for the nucleus ^°Ca for the various orbital 
angular momenta/ in an energy interval [-50 MeV, -10 MeV]. The experimental data 
are from ref. [17]. Note the different scales for the spectral function in the various 
parts of the figure 
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